Spontaneous solitons in the thermal equilibrium of a quasi-one-dimensional Bose gas 
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We show that solitons occur generically in the thermal equilibrium state of a weakly-interacting elongated 
Bose gas, without the need for external forcing or perturbations. This reveals a major new quality to the experi- 
mentally widespread quasicondensate state, usually thought of as primarily phase-fluctuating. Thermal solitons 
are seen in uniform ID, trapped ID, and elongated 3D gases, appearing as shallow solitons at low quasiconden- 
sate temperatures, becoming widespread and deep as temperature rises. This behaviour can be understood via 
thermal occupation of the Type II excitations in the Lieb-Liniger model of a uniform ID gas. Furthermore, we 
find that the quasicondensate phase includes very appreciable density fluctuations, while leaving phase fluctua- 
tions largely unaltered from the standard picture derived from a density-fluctuation-free treatment. 



Solitons, or non-destructible local disturbances, are im- 
portant features of many one-dimensional (ID) nonlinear 
wave phenomena. In ultra-cold gases, they have long been 
sought, and were first observed to be generated by phase- 
imprinting HI More recently, their spontaneous formation 
in ID gases was predicted as a result of the Kibble-Zurek 
mechanism 0], rapid evaporative cooling J^t], and dynam- 
ical processes after a quantum quench^. Here we show 
that they actually occur generically in the thermal equilibrium 
state of a weakly-interacting elongated Bose gas, without the 
need for external forcing or perturbations. This reveals a ma- 
jor new quality to the experimentally widespread quasicon- 
densate state. It can be understood via thermal occupation of 
the famous and somewhat elusive Type II excitations in the 
Lieb-Liniger model of a uniform ID gas H. 

A mathematically distinct class of soliton equations are 
the completely integrable systems. Among them, the Gross- 
Pitaevskii[8, 9] equation describes weakly interacting bosons 
in a ID geometry in the mean field approximation. The cor- 
responding multi-atom Lieb-Liniger model of N bosons on 
the circumference of a circle interacting by contact forces Jioll 
has elementary excitations of two kinds: those of a Bogoli- 
ubov type and an additional "Type II" branch [7]. These ad- 
ditional excitations have been associated with solitons of the 



mean field II 1114 1411 . Although the trapping potential removes 
integrability, from the early days experimenters have searched 
for gray solitons. See lfl5ll for a review. 

Typically, by irradiating one part of the condensate, one en- 
gineers a phase difference with the remainder, and a dark soli- 
ton forms at the interface between the phase domains JH 
Other proposed schemes involve taking the system away from 
equilibrium J31-0]. Our results show that solitons are in fact 
present spontaneously even in equilibrium. However, engi- 
neered solitons have been easier to identify with the standard 
destructive imaging measurements because their position does 
not vary from shot to shot. 

Here we generate a classical field ensemble that describes 
the weakly interacting Bose gas at thermal equilibrium 111611 . 



and then show that gray solitons indeed are already there. This 
is demonstrated by tracing the time evolution of a single copy 
of the system. We find spectral properties consistent with Lieb 
Type II excitations. Moreover, we will also show that the pres- 
ence of gray solitons in thermal equilibrium remains valid for 
very elongated traps that are no longer strictly ID. 

Firstly, let us consider the simplest situation: a uniform ID 
weakly-repulsively-interacting gas in free space with periodic 
boundary conditions, as per Lieb and Liniger To obtain 
the equilibrium state, we apply the classical field approxima- 
tion (CFA) Inl [TH] . Within this approach the usual bosonic 
field operator ^f(z) which annihilates an atom at point z is 
replaced by an ensemble of complex wave functions ^Sf(z) ob- 
tained using a Monte Carlo sampling algorithm |[l6|] (see Sup- 
plementary material for details). 

In Fig. Q] we show the time evolution of the density and 
phase for single realizations in equilibrium at several temper- 
atures. The starting conditions are randomly chosen from the 
collection of wavefunctions for the canonical ensemble ob- 
tained in the CFA. Subsequent time evolution of these chosen 
realizations in the CFA is via the Gross-Pitaevskii equation of 
motion. 

The only relevant dimensionless parameter appearing in 
the Lieb-Liniger model of a ID interacting Bose gas is 7 = 
mg/h 2 p, where g characterizes the atom-atom interaction 
and p = N/ L is the linear density of the system. In what fol- 
lows we will use L, mL 2 /h, H 2 /mL 2 , and H 2 /mL 2 kB as the 
units of length, time, energy, and temperature, respectively. 

Firstly, the top panels show the system at the relatively low 
T = 10 4 = /i/2. There are numerous density disturbances in 
the form of dips and peaks, traveling in both directions, near 
the speed of sound (which is depicted by the red line). There 
exist both defects that travel faster, and those that travel slower 
than sound. The deepest dips are slower - the slope of their 
trajectories on the plot is higher than that of the red line. At 
low temperature only fast moving shallow defects are present. 
Some of them correspond to packets of Bogoliubov phonons, 
others are solitons. This can be checked by inspection of the 
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the accepted dips is compared with the Lieb Type II excita- 
tion spectrum in Fig. [2] for three temperatures. At the low- 
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FIG. 1 . Density (left) and phase (right) as a function of time of a 
ID Bose gas at equilibrium for a single realization. Here 7 = 0.02 
(N = 10 3 ) and the temperatures are T — 10 4 = /jl/2 (top), 7 x 10 4 
(middle), and 10 5 (bottom) in units of h 2 /mL 2 kB- The red line in 
the upper left panel corresponds to travel at the speed of sound. 



phase jump across the defect (for a soliton its sign is related 
to its direction of movement), and by fitting the local density 
dip to a soliton solution [19]. However, shallow fast solitons 
can be hard to tell from phonons. 

For a higher temperature T = 7 x 10 4 (middle panels), 
long-lived deep dips travelling far slower than the speed of 
sound appear. These become prolific as the temperature in- 
creases to T = 10 x 10 4 (lower panel). An analysis of the pa- 
rameters of slower- than- sound defects confirms their interpre- 
tation as dark solitons. Indeed, one sees the anti-correlation 
between the depths and speeds of the dips that is expected 
for dark solitons. Near maximal depth the dips approach be- 
ing stationary. One also sees numerous soliton collisions, and 
associated phase-shifts, with the number of solitons usually 
conserved. 

The match between dark solitons and Lieb Type II exci- 
tations is seen by comparing fitted parameters describing the 
local density dips to the grey soliton solutions. To do this, 
we fit a density profile (characterised by soliton depth, local 
density and soliton position) to each dip. If the rms parame- 
ter errors are smaller than 10% of the local density (for depth 
or density) or the healing length (for position) we accept the 
dip as a soliton. We then calculate soliton energy e and mo- 
mentum p from the parameters lUlll . collecting data from many 
time snapshots. The resulting energy-momentum relation of 




FIG. 2. Comparison of the spectra of Lieb Type II excitations (solid 
black line) and dips in the thermal state for three different tempera- 
tures: T = 10 4 (red dots), 4 x 10 4 (green dots), and 10 5 (blue dots). 
Here 7 0.02 (N 10 3 ). 

est temperature, the spectral match is good, and includes the 
negative curvature of the type II branch. For higher tempera- 
tures, the fitted soliton energy lies somewhat below the Lieb 
II line, retaining a qualitative match. Interpretation of the shift 
is difficult as interactions between solitons are frequent at this 
temperature. Notably, one sees a filling of the high p spectrum 
with temperature. 

This allows us to explain the temperature dependence of the 
appearance of deep solitons. These have the maximum Lieb II 
energy and momentum Jl3ll . For the parameters of Fig. [TJ that 
is about 10 5 in units of temperature (see Supplementary ma- 
terial). Thus we would expect the deepest solitons to appear 
when temperatures of this order are reached. This is what is 
seen. We found similar agreement for 7 = 0.002. 

It is appealing to think of the system as of a gas of quasipar- 
ticles of two kinds - where two bosonic excitation families co- 
exist like in the Lieb-Liniger model (see Supplementary ma- 
terial). These two kinds - phonons and solitons, interact with 
each other in an inelastic way, so that solitons can be born 
from Bogoliubov excitations and vice versa. Such interactions 
could allow the gas of quasiparticles to reach the equilibrium 
state. 

Now a question that begs to be asked is whether solitons 
will also be present in the less idealized case of a trapped gas 
in the quasicondensate regime - a situation ubiquitous in con- 
temporary experiments. The CFA can be applied in a simi- 
lar way as above (see Supplementary material). Firstly, con- 
sider a harmonically-trapped ID gas. Fig. [3] (left) shows its 
density evolution for two different temperatures. The lower 
temperature case (left upper panel) has T chosen such that 
the phase coherence length 1$ is approximately equal to the 
size of the cloud W, i.e. it corresponds to the temperature 
which separates the quasicondensate and true condensate 
phases, as first described by Petrov et al. [20]. For temper- 
atures lower than this, only shallow fast solitons are seen in 
the bulk of the cloud. At much higher temperature (lower left 
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panel in Fig. [3]) a mass of deeper long-lived solitons emerges. 
They remain within the main part of the atomic cloud, and 
their number and depth increases with temperature. One can 
easily check that for the lower temperatures the solitons oscil- 
late in the trap with a frequency which within a few percent 
agrees with uj / \[2 from a simple model Q2JJ] • At high tempera- 
tures, inter-soliton collisions come to dominate the oscillating 
behavior, making that model inappropriate. 
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FIG. 3. Solitons in trapped clouds. The left panels show the time- 
evolution of a trapped ID condensate with parameters: N = 1000, 
ID interaction strength g = 0.31 in units of huo z (h/m(jO z ) 1 ^ 2 , 
and trap frequency uj z — 2tt x 10 Hz. The temperatures are: 
}zbT — lbhwz (upper frame) and fcsT — 260/^ (lower frame). 
The right panel shows the evolution of the central density in a three- 
dimensional elongated trapped cloud (horizontal direction). Here, 
the trap frequencies are: uj z = 2ty x 10 Hz and uj± = 2tv x 1000 Hz, 
with again N = 1000, and the (3D) interaction strength g^D — 
0.213 in units of fluj z (h/m(jj z ) 3 ^ 2 . The temperature is 80 Huj z /kB- 



The center of the cloud can be considered as a uniform gas 
under a local density approximation (LDA) as confirmed in 
recent experiments 11221 12311 . This is described by the dimen- 
sionless interaction strength (7 = g/(n}) and temperature 
(r = 2k B Tm/h 2 (n) 2 ). In Figures and ffl the T = 15 
gas lies deep in the quasicondensate regime (7 = 0.0033, r = 
0.0034), while the center of the T = 260 cloud is a decoherent 
quantum gas Hi (7 = 0.0059, r = 0.19). Here, FigHshows 
that significant local density fluctuations Sn/ (n) « 10 — 40% 
are present in the quasicondensate phase r <C y/j. This is 
to be contrasted with the common view that the only notable 
fluctuations in the quasicondensate are those of the phase, a 
situation we see only at the very lowest temperatures. (Note 
also the long- wavelength density fluctuations in ll23lo . 

Importantly, the CFA fluctuations agree with predictions 
obtained from the exact local density correlation function 11241] 
in the Yang & Yang description! 25] of the interacting uni- 



form gas. Note that we have taken account of the fact that 
the CFA does not include zero temperature fluctuations, so 
that the Sn in Fig. 01 shows only the thermal contribution. 
The equivalent contribution in the Yang and Yang solution is 

(Sn) therms /(n) = y/ gM (0) - #T 2 io( )' wmcn is shown as 
a black line. The level of agreement implies that the solitons 
we have seen here in Fig. [3 that contribute a large proportion 
of the density fluctuations, must also be an essential feature of 
the Yang&Yang description. 




FIG. 4. Phase and thermal density fluctuations in the center of the 
trapped ID gas as a function of temperature. Parameters are as in 
Fig. [3] except for varying T (in units of Huj z /kB). Purple circles 
show Sn/ (n) in the center of the trap; Orange squares the ratio of the 
phase coherence length 1$ to the cloud width W. The latter is taken 
to be the distance between points at which the density has fallen to 
10% of the central value. Grey bars indicate temperatures T</> and an 
estimate of the location of the quasicondensate / decoherent quantum 
gas crossover at r = y/j 1 24] . Yang&Yang predictions (black solid 
line) are for the average density in the center (\z\ < 2) of the cloud. 
Some statistical uncertainty is visible. 

We have also calculated the coherence length 
through matching the phase correlation function 
(z, z') = p\ (z, z') I \J ' pi(z, z)pi(z' : z') to an exponential 
decay exp [— \z — z'\/l^\. This is shown in Fig. [4] in orange, 
and is still in agreement with the canonical treatment that 
ignores density fluctuations lEoll . The one-particle density ma- 
trix, is obtained by an averaging of (z)^(z')) 
over the initial canonical ensemble. For the low temperature 
case (T = 15), 1$ = 13.7a ho = 13.7 y/h/u z m, while 
the cloud width W is 14.7aho- To compare, the standard 
quasicondensate expressions in lEoll that consider only phase 
fluctuations give an estimate of the temperature at which 
1+ = W. It is TP haseonly = 3N 8 ( ^f = 12.5 huj z /k B for 
these parameters, which is a good agreement with the CFA 
calculation. Our agreement on with the standard descrip- 
tion is so good because it is a manifestation of the fact that 
the correlation length obtained from a pure phase dependence 
of the classical field, i. e. g^\z, z') = exp[z(0(z) — (j){z'))\ 
like in 112011 — where (j)(z) is the phase at point z — does not 
differ significantly from what is obtained via the complete 
correlation function g( l \z,z') [26] . We conclude then, 
that the 10-40% density fluctuations and solitons that we 
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see in the quasicondensate, are not inconsistent with the 
past calculations of phase fluctuations. Rather, they are an 
inherent feature that had remained un-noticed until this time 
due to the prevalent focus on phase coherence. 

These results raise an interesting issue regarding the 
Kibble-Zurek mechanism of defect formation (KZM) in ID 
gases. In a recent paper H it was conjectured that dark soli- 
tons should be created by the KZM during rapid cooling of the 
gas(27l|28]. Simulations Q] of the cooling of a ID Bose gas 
in a harmonic trap indeed revealed the presence of solitons. 
However, we know now that some solitons are present even 
in thermal equilibrium. Only in the case when the number of 
solitons at the end of the cooling process exceeds that found 
in its corresponding thermal equilibrium, could simulations of 
that kind confirm the KZM at work. 

Finally, there is a question whether these solitons can still 
survive when the strict ID trapped system crosses over to one 
in a very elongated cigar-shaped trap. This is made plausible 
by past work which showed that quasicondensate-BEC-like 
transitions occur for cigar-shaped systems analogously with 
the true ID gas!29l-[311. The right panel of Fig. [3] shows the 
evolution of a thermal state in a fully three-dimensional cal- 
culation. Solitons are still clearly visible. This indicates that 
they may be present, and indeed even widespread, in many 
existing experiments. 

The observation of thermal solitons with existing equip- 
ment has been non trivial because of two factors: The soli- 
ton width in- situ was typically significantly narrower than the 
detector resolution, and an identification of solitons by "eye- 
balling" is only straightforward if one has access to observa- 
tions of the in-situ dynamics, not destructive snapshots. Look- 
ing at the edges of the gas offers an indirect way of detecting 
the solitons B32I1 . As they bounce between the edges of the gas, 
they approximately retain their absolute depth while also in- 
creasing in width and reaching close to 100% relative depth at 
the turning point, making them accessible with realistic reso- 
lutions. We do not see correspondingly deep phonon distur- 
bances in our images. Therefore 113 211 . we have counted the 
density of dips with over 90% relative depth in the wings of 
the ID cloud. The growth of the density of these solitons with 
temperature is shown in Fig. \5\ 

Further progress could be made by recent experimental ad- 
vances such as high-resolution detection based on scattering 
of an electron beam l|33[l , analysis of momentum spectrafl, 
phase-jump statistics [Ml (phase domains are larger than soli- 
ton dips) or with long Ol 0, 35] or anti-trapped expansions fl36ll 
(See Supplementary material for more detail). 

In summary we have found that solitons are a natural 
and spontaneous feature of the ID and elongated weakly- 
interacting Bose gases. While solitons have been studied in 
such systems before by imprinting, the above results show 
that they are in fact common. The good agreement with ex- 
act density fluctuations in Fig. [4] shows that quantitative study 
of these phenomena is possible with the CFA, and possibly 
related methods. The ID Bose gas at thermal equilibrium is 
a system where two very different kinds of bosonic excita- 
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FIG. 5. Number of solitons per /jm on the edges of the cloud, per 
snapshot, versus temperature. There are 10 4 Rb 87 atoms in a ID 
trap with oj z = 2tv x 1.9Hz. The ID coupling assumes a Gaussian 
transverse profile corresponding to lo± = 2tt x 128Hz. Counting 
was over 7 /mi intervals on either side of the cloud centered around 
the position where the average density is 10/xm -1 (12% of peak). 
Dips whose central density falls to below 10% of the mean density 
at their trap position were counted. 



tions are present simultaneously. The system actually con- 
tains much more than the standard picture, especially in the 
quasicondensate regime. Here, in addition to the well known 
phase fluctuations, appreciable density fluctuations are found, 
including deep solitons. These spontaneous solitons in the 
thermal state are analogous to pairs of vortices present in a 2D 
gas near the Berezinskii-Kosterlitz-Thouless transition l37ll . 
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Supplementary material for: 
Spontaneous solitons in the thermal equilibrium of a quasi-one-dimensional Bose gas 



The classical field approximation 

Within this approach HI 0] the usual bosonic field opera- 
tor $(z) which annihilates an atom at point z is replaced by 
the complex wave function S&(z). Technically speaking, we 
first expand the field operator in the basis of one-particle wave 
functions, appropriate for the problem considered. Then, ex- 
tending the original Bogoliubov idea [3] to all macroscopi- 
cally occupied one-particle modes, we replace the operators 
corresponding to these modes by c-numbers. Restricting the 
expansion only to these modes, the field operator is turned into 
a complex wave function - the classical field. For a plane- 
wave basis, 



*(*) 



\k\<k n 



(1) 



The above expression is appropriate for bosons confined in a 
box of length L with periodic boundary conditions. For a har- 
monic trapping potential, the single-particle trap eigenfunc- 
tions are more convenient. The summation is extended over 
all modes up to the momentum cut-off hk max . The optimal 
choice of the cut-off is discussed in Ref . [0] . Fully three di- 
mensional, elongated trap simulations require recalculation of 
the optimal cut-off condition since the explicit results in [0] 
are valid only for symmetric D-dimensional traps and, more- 
over, are asymptotic for large number of atoms. Each classical 
field, ^f(z), shares many properties with the single-shot mea- 
surements of the atomic cloud that occur in experiment. The 
classical field satisfies the following equation of motion 0JLD : 

zft^*CM)= (-^^ + gmz,t)\ 2 )9(z,t), (2) 



dt 



where g characterizes the atom-atom interaction and the 
nonlinear term is projected on the subspace spanned by the 
macroscopically populated modes. 



Thermal states 

To obtain the thermal equilibrium state of a ID Bose gas 
within the CFA, we numerically generate members of the 
canonical ensemble of states [5], i.e. states populated accord- 
ing to the probability distribution given by: 

P({a k }) = L- E */ k ° T , (3) 



where Z is the canonical partition function and T is a temper- 
ature. The energy, Ey, accumulated in the classical field is 
given by : 

E* = J dz**(z)H *(z) + \d J l dz\^{z)\ A . (4) 
where Ho is the single-particle Hamiltonian 

h 2 d 2 



Hq = T4 r ap(z) 



2m dz 2 



(5) 



with the trapping potential V trap (z). An extra constraint on 
the amplitudes should be fulfilled: 



(6) 



\k\<k n 



where TV is the number of atoms. An efficient way to ob- 
tain states belonging to the canonical ensemble at given tem- 
perature T is a Monte Carlo method using the Metropolis 
algorithm[6]. Here, a random walk in the phase space of the 
system is performed and all visited states become the mem- 
bers of the canonical ensemble. These states are used to cal- 
culate statistical averages of any observable. Details for ultra- 
cold Bose gas systems are given in Ref. I5fl. 



Temperature at which deep solitons appear 

Assuming the excitations observed within the CFA are Lieb 
type II excitations, we are able to estimate the temperature at 
which deepest solitons appear in the thermal state of a uniform 
system. To obtain the dispersion curve for the type II excita- 
tions, one is required to solve the inhomogeneous Fredholm 
integral equations J7D . This can be done numerically, and for 
7 = 0.02 the dispersion curve, i.e. energy versus momentum 
e(p), is found to be well approximated by e(p) = app + bp 2 , 
where a = 0.125 and b = —0.021. Maximal excitation energy 
occurs for the highest momentum on the curve p max = ttN 
(see [lZD for details), and corresponds to the deepest, stationary 
solitons. It is then given by the formula e max = 7rN 2 (a + irb) 
in our units, which for N = 1000 is about 10 5 . Therefore, for 
7 = 0.02, the deepest solitons are expected at temperatures of 
this order or higher. 
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Spectral analysis 

Further arguments for the simultaneous presence of both 
kinds of excitations (Type I - Bogoliubov phonons, and Type 
II - solitons) in a ID Bose gas at the equilibrium come from a 
spectral analysis of the classical field. Fig. [T] shows the space- 
time spectral density \^f(u), k) | 2 of a single realization at equi- 
librium. Like in the three-dimensional case (see Ref. [8]), two 
curves crossing at a frequency equal to the chemical potential 
ji = 2 x 10 4 are clearly visible. The low-momentum slope of 
each curve is just the speed of sound. The parts of these curves 
below the chemical potential are necessary to construct the 
phonon-like section of the excitation spectrum [8,9] and fade 
at momenta at which the dispersion curve changes its charac- 
ter from linear to parabolic. This part of the spectrum proves 
the existence of Bogoliubov phonons in thermal equilibrium. 
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Figure S 1. Spectral density of a ID Bose gas at equilibrium for 
7 = 0.02, T = 10 4 within the classical field approximation (up- 
per panel, detail - lower left panel). In addition to the Bogoliubov 
phonon nature of the main curves, the spectral density exhibits a 
checkerboard pattern for low momenta. This pattern is related to the 
presence of solitons in the system. Bottom right panel: spectral den- 
sity of two dark solitons propagating in opposite directions. Soliton 
velocities are ±0.8c, with c the speed of sound. The checkerboard 
pattern again appears on top of the main structure consisting of two 
lines uj = ±uk. It is interpreted as a result of interference between 
solitons. 



board pattern appears as a result of interference between the 
solitons, and thus it is a strong signature of the existence of 
dark solitons in the system. Indeed, the checkerboard pattern 
becomes more regular when the number of solitons increases. 



The classical field approximation in one dimension 

The classical field approach|[l|] and the very closely related 
PGPE (Projected Gross-Pitaevskii Equation) approach^ 
have been benchmarked on numerous occasions — see de- 
scription in the reviews JH IH, or lUlll . To demonstrate its 
applicability to the density fluctuations of ID gases as consid- 
ered in this article, we have compared CFA results to a recent 
experiment that measured density fluctuations [1121. 

Using the actual parameters of Jl2ll and the prescription for 
the cut-off parameter used in this article, we obtained the very 
good agreement illustrated in Fig S|2l 
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Figure S 2. Local atom-number variances, S 2 n, as a function of 
the mean local atom density in a weakly interacting quasi ID gas 
confined by a harmonic potential. Red circles are taken from Fig. 
lc. in I12n . Blue crosses are the results of our ID simulation 
with parameters from the experiment IU2I1 : g±D — 2huu±a = 
2.06huo z (h/m(jj z ) 1 ^ 2 is the coupling constant, a = 5.7nm is the 
3D s-wave scattering length, and u± = 2tt x 3.9 kHz, u z = 2tt x 4 
Hz are the frequencies of the transverse and longitudinal harmonic 
confining potentials, respectively. T = 0.09ftw± /kB = 88fiw z /kB. 



However, in addition to the phonon-like behavior, the 
spectral density exhibits a phenomenon absent in the three- 
dimensional system - in the region of low energies and mo- 
menta a checkerboard pattern appears (see Fig. Q]). This is a 
signature of the presence of solitons. We confirm this inter- 
pretation by considering the spectral density of two counter- 
propagating dark solitons moving with the same velocity u 
(and depth), shown in the lower right panel of Fi g. ED These 
solitons are obtained from the Zakharov solution 111 PL and its 
subsequent evolution according to the nonlinear Schrodinger 
equation. Finally, the dispersion curve of a single dark soliton 
consists of only a single line uo = uk. Hence, the checker- 



Possible detection schemes 

Let us now turn to the question of experimental observation 
of solitons in the thermal state. Unfortunately, they are not 
straightforward to directly detect because of two factors: (1) 
The soliton width in-situ, being a fraction of a /im, is signif- 
icantly narrower than the usual detector resolution (typically 
several /im). However, a recently developed detection scheme 
based on scattering of an electron beam could overcome this 
difficulty Jl3ll . Also, (2) an identification of solitons by "eye- 
balling" is only straightforward if one has access to observa- 
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tions of the in-situ dynamics - a single-time density slice typi- 
cally shows many density dips, and it is not easy to distinguish 
the fast "phonon-like" ones from bona-fide solitons. 
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Figure S 3. Experimental signatures of thermal solitons in the uni- 
form gas — phase-jump frequency (log scale) at 7 = 0.02 for low 
temperature (no deep solitons - blue, T — 10 4 ) and high tempera- 
ture (many deep solitons - red, T — 10 x 10 4 ). Columns show the 
effect of different imaging resolutions in units of healing length £ ( 
= 1/Ny/j in our units), which is the half- width of a deep soliton. 
Resolution worsens from left (below £) to right (unable to resolve 
soliton dips). Note the strong resistance of the difference in shape 
of the distribution (parabolic / extended) even into regimes where 
individual solitons cannot be resolved. 

However, several approaches hold promise of overcoming 
one or both of the above issues. One is a direct observation of 
deep solitons in an expanding cloud. After about 10ms of ex- 
pansion of a typical cloud, the soliton size can exceed detector 
resolution. The method was successfully used in engineered 
soliton experiments 

Another promising approach is to look for the expected 
large phase jump between phase domains that occurs at the 
soliton. If relatively recognizable phase domains are present 
between solitons, then the phase jump of close to ±tt that 
occurs at deep soliton defects should be detectable with an 
imaging resolution that is sufficient only to resolve the phase 
domains 111711 . These are much wider than the solitons if the 
soliton density is not extreme. 

It turns out that even with very noisy domains, a qualitative 
difference can be observed between clouds with and without 
solitons. Fig. S(3] shows the phase-jump histogram (for phase 
jumps between neighboring pixel pairs) at 7 = 0.02 for both 
low temperature and high temperature. The low- temperature 



system has a Gaussian distribution of phase jumps (parabolic 
on the plot, which has a log scale), which broadens but does 
not change shape as the resolution is worsened. This reflects 
the addition of more and more random small phase fluctu- 
ations that arise primarily from Bogoliubov excitations. In 
stark contrast, when deep solitons are present, phase jumps of 
7r appear quite frequently, and - most significantly - the distri- 
bution flattens out for large jumps. This qualitative difference 
(flattening-out/ Gaussian) survives even to resolutions that are 
incapable of resolving the actual density dip of the solitons. 

The anti-trapping technique may also be quite promising. 
This technique has been used to image short wavelength 
shock waves in elongated gases fl8ll . which have some phe- 
nomenological features in common with solitons. An inverted 
parabolic potential is rapidly applied to the gas along the long 
direction to prevent the evolution of in-situ density to momen- 
tum density in that direction, but to instead magnify the den- 
sity profile. 

Finally, a recent work has presented signatures of soli- 
tons in the form of characteristic features of the momentum 
distribution Jl9ll . It analyzed solitons formed after a quantum 
quench, but is also appropriate for the randomly placed ther- 
mal solitons discussed here. 
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